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ABSTRACT 

Semi-analytic models of galaxy formation postulate the existence of virialized 
gaseous halos around galaxies at high redshifts. A small fraction of the light emitted 
by any high-redshift quasar is therefore expected to scatter off the free electrons in the 
halo of its host galaxy. The broadening of the scattered emission lines of the quasar can 
be used to measure the temperature of these electrons. For gas in virial equilibrium, the 
velocity width of the scattered line profile is larger by a factor of ~ {mp/meY^'^ — 43 
than the velocity dispersion of the host galaxy and reaches ^ 10, 000 km s~^ for the 
massive galaxies and groups in which bright quasars reside. In these systems, the 
scattered width exceeds the intrinsic width of the quasar lines and hence can be used 
to measure the virial temperature of the quasar host. The high degree of polarization 
of the scattered radiation can help filter out the extended scattered light from the 
central emission by the quasar and its host galaxy. The signal-to-noise ratio of the 
spectral broadening can be enhanced by matching the full spectrum of the scattered 
radiation to a template of the unscattered quasar spectrum. Although the central fuzz 
around low-redshift quasars is dominated by starlight, the fuzz around high-redshift 
quasars might be dominated by scattering before galaxies have converted most of their 
gas reservoirs into stars. 

Subject headings: galaxies: fundamental parameters - quasars: emission lines 



1. Introduction 

The depth of galactic potential wells has been traditionally inferred by dynamical methods 
(HI or stellar rotation curves and satellite dynamics) or gravitational lensing. Unfortunately, the 
application of these methods to galaxies at very high redshifts proves difficult. The resulting 
uncertainty in the velocity dispersion of galactic halos prevents current data sets on the abundance 
and clustering of galaxies at z ~ 3 from being able to discriminate among different cosmological 
models (Adelberger et al. 1998). 

Semi-analytic models of galaxy formation and evolution (e.g., Kauffmann, White, & 
Guiderdoni 1993; Baugh et al. 1998; Somerville, Primack, & Faber 1998) postulate the existence of 
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gaseous halos at the virial temperature around galaxies at high redshifts. The halos are expected 
to contain gas that was assembled into the galaxies but was unable to cool between successive 
heating episodes, caused by either merger events or supernova-driven winds. The supply of cold 
gas from the dense core of these halos controls the star formation rate of their host galaxies. 
Although X-ray halos are rare among local galaxies (Bregman & Glassgold 1982; Fabbiano 1989), 
they might have been prominent during their early evolution. A direct measurement of the gas 
temperature in these halos could determine the potential depth of their host galaxy. 

In this Letter we propose a novel method for probing the virial temperature of galactic halos, 
in those cases where a bright quasar shines at their center. The quasar radiation is scattered 
by the free electrons in the halo and the scattered emission lines are Doppler broadened by the 
thermal velocity of these electrons]^. Due to the small electron mass, the velocity width of the 
scattered lines is larger by a factor of ~ (mp/mg)^/^ = 43 than the velocity dispersion of the 
galaxy, and could reach values ^ 10, 000 km s~^ for the massive galaxies and groups in which 
bright quasars reside. In these objects, the scattered line width exceeds the typical intrinsic width 
of the quasar emission lines (FWIIM~ 5,000 km s~^; see, e.g. Peterson 1997 for a review) and can 
be used to measure the electron temperature. Since quasars tend to reside in dense environments 
such as groups of galaxies (Hartwick & Schade 1990; Bahcall & Chokshi 1991; Ellingson, Green & 
Yee 1991; Fisher et al. 1996), additional scattering may also result from the intra-group gas. 

In §2 we calculate the expected shape and amplitude of the scattered quasar lines as a 
function of the gas temperature and density profile. Finally, §3 discusses the implications of our 
results. 



2. Scattered Line Profile 

The redistribution function n'; z^, n)di^'(iz^(dil'/47r)((i$7/47r) provides the probability of 
scattering a photon from a frequency between {v', v' + dv') and solid angle (Kl' in the direction n' 
to a frequency between [v, u + dv) and solid angle dil in the direction n. For Thomson scattering 
by a Maxwellian population of electrons at a temperature T, the redistribution function is given 
by (Dirac 1925; Mihalas 1978), 
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where = n' • n = cos 9, and 



2kT \ 
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^Interestingly, this broadening effect was only recently detected in a local astrophysical setting - the solar wind, 
by the SOHO satellite (Fineschi et al. 1998). 
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is the typical thermal velocity of an electron in units of the speed of light c. Here rrie is the 
electron mass, and k is Boltzmann's constant. The redistribution function can be integrated over 
any arbitrary scattering geometry to obtain the profile of a scattered quasar line. 

As an illustrative example, we consider below the case of a spherical, optically-thin, gas cloud 
which scatters a narrow line of frequency z^o ^ind total luminosity Lq from a central source. We 
assume that the source luminosity per unit frequency is L^i = Lq6{u' — uq). An electron at a 
radius r scatters a net fraction ((TT/47rr^) of Lq, where cjt = 0.67 x 10^"^^ cm^ is the Thomson 
cross-section. Hence, the scattered spectral luminosity at a frequency v is 

L^ = J ne (^) l^l'^Jd^' R{^',n';u,n)LA^'), (3) 

where ne(r) is the number density of electrons at a radius r. By substituting d^r = Airr'^dr and 
Ljy/ = Lq6{v' — uq) and integrating by parts, we get (Rybicki &: Hummer 1994) 

-A'^IA /q a2 a4\ IAI /IAI\1 

(4) 
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where, erfc(x) = {2/y/T:) /^exp(— is the complementary error function; 

A.<^. (5) 

is the frequency shift relative to the line center in units of the typical Doppler shift of a thermal 
electron, and 

(•OO 

dr UeaT, (6) 



is the radial optical depth of the halo. The integral of the right-hand-side of equation (^) over A 
is unity. The scattered profile of a line with a finite initial width can be obtained by convolving 
the initial line profile with the kernel given by equation (Q). This kernel is plotted as the solid line 
in Figure 1. 

The effective width of the scattered line can be found from the second- moment, 
(A2) = {J A'^UdA)/ J L^A. We find. 



V \ vo 

where the angular brackets denote an average over L^. 

Next, we apply the above results to the particular halo profile of an isothermal sphere 
with a one-dimensional velocity dispersion a^. For a fully-ionized gas of electrons and protons, 
kT = rupal, where nip is the proton mass. We then get from equation (^), 

-^c = 1.2 X 10^(7200 km s"\ (8) 



- 4 - 



where o"2oo = (o"^/200 km s~^). We assume that the gas provides a fraction /gas of the total mass 
and follows the radial profile of the dark matter at large radii, but has a core of radius re- Hence, 
we take = /gasC^/[2vrGmp(r^ + r^)], where G is Newton's constant. Equation (|6|) then yields 
the net fraction of the quasar light that gets scattered, 

r = 2 X W~^hoaloon\ (9) 

where /lo = (/gas/10%) and ri = [vc/l kpc). The value of /gas can be as high as 20% if the 
cosmological density parameters of baryons and dark matter are 0.06 and 0.3, respectively. The 
radius inside which the cooling time due to free-free emission is shorter than a fraction rj of the 
Hubble time is ~ 25 kpc (?7/ioo"20o)"'^^^[(l + z)/^~^^^. The relevant numerical value of ^ 1 
depends on the average fraction of the Hubble time that passes between major mergers, as 
these mergers could reheat the diffuse gas to its virial temperature. In fact, if quasar activity 
is triggered by galaxy mergers, then r] may be as small as the ratio between the dynamical 
timescale of the host galaxy (during which the quasar is fueled) and the Hubble time, namely 
~ 0.1(r/30 kpc)(T2oo[(l + z)/A]^/'^. Inside the cooling radius, hot gas may still be replenished by 
supernova-driven winds (in which case the geometry of the scattering material is likely to be 
non-spherical). The quasar radiation itself is expected to ionize and heat much of the quiescent 
interstellar medium as well, and so overall it is natural to expect ri ^ 10. In resolution-limited 
observations where the quasar has to be subtracted from the image, the effective value of r would 
be dictated by the extent of the quasar mask which sets a minimum value to r. Scattering by dust 
would also contribute to the diffuse light, but would not broaden the quasar lines significantly. For 
the massive galaxies and groups in which bright quasars tend to reside (Hartwick & Schade 1990; 
Bahcall & Chokshi 1991; Ellingson et al. 1991; Fisher et al. 1996), the width of the scattered 
emission lines will typically be ^ 10^ km s~^, i.e. in excess of their intrinsic width. 

Note that the optical depth of a cluster of galaxies with ~ 10^ km s~^ and a core radius 
Tc ~ 250 kpc is the same as that of a galaxy with ~ 200 km s~^ and Vc ~ 10 kpc. However, 
in the cluster case, the surface brightness of the scattered light is a factor of 25^^ = 625 smaller 
and is much more difficult to separate from contaminating background light. Moreover, the lines 
scattered by the cluster would only appear as a broad hump on top of the continuum because 
of their large width, (bv jv) ^ 0.2. Scattering of beamed quasar light by intracluster gas was 
suggested as the origin of the observed alignment between the optical continuum and radio axes 
of high redshift radio galaxies (Fabian 1989). 

The signal-to-noise ratio for the detection of electron broadening can be improved by using 
a full template of the quasar spectrum with all of its lines, such as O VI, Lya, C IV, Mg II, and 
Ha, /3, 7. The scattered lines will be centered on the template wavelengths and be broadened by 
the same amount (except for the blue wing of the Lya line which gets absorbed by the Lya forest 
of the intergalactic medium). The broadening due to electron scattering can be distinguished from 
the intrinsic width of the lines since it should be the same for all lines while the intrinsic width is 
often different for different lines (Peterson 1997). 



- 5 - 



For a density profile rig cx r~° the degree of polarization near the limb of the scattering halo 
is P = (a + l)/(a + 3) (Cassinelli & Hummer 1971). In the singular isothermal sphere case, a = 2, 
and so P = 0.6. Since the intrinsic quasar emission is often weakly polarized, this high degree of 
polarization can be used to filter out the scattered component from the emission component in 
the spectrum of the quasar and its host (and may also be used to reduce the surface brightness 
contrast between the quasar pointlike emission and the extended scattered light). Since the net 
polarization of a spherical halo is zero, the detection of polarization requires that the halo be 
resolved. At any point, the E-vector polarization axis is perpendicular to the radius vector from 
the halo center. 

In spectroscopic observations for which the halo is resolved, the line profile might deviate 
from the shape derived above. In particular, the spectral shape in equation (^) holds only if the 
flux is summed over a symmetric ring or over half of the symmetric projection of the entire halo on 
the sky (with the value of r depending on geometry). If a spectrum is taken for any other portion 
of the halo, then the line profile would be different. For example, in the limit of high angular 
resolution (needed for the detection of polarization) and a steep electron density profile (for which 
the polarization approaches 100%), the line shape would approach R{vq,i') for a scattering angle 
oi = 90° (/X = 0), since most of the scattered flux would be contributed from the point of closest 
approach of the line-of-sight relative to the source. The dashed line in Figure 1 compares the 90° 
scattering profile of equation (||) to the spectral shape given by equation 



3. Discussion 

We have shown that the scattered light of a central quasar can serve as a thermometer of 
hot galactic halos with a virial temperature T = 5 x lO^cloo K. This probe complements other 
methods for measuring the depth of galactic potential wells, although it is limited only to galaxies 
which harbor a quasar]^. Local observations of galactic centers (Magorrian et al. 1998) and popular 
models of quasar evolution (e.g., Haehnelt & Rees 1993; Haiman & Loeb 1998) suggest that all 
galaxies went through a quasar phase during their history, and so this method could be applied to 
a fair sample of all galaxies. 

The scattered flux is expected to replicate the quasar spectrum with broadened line widths of 
~ 10^(T2oo km s~^, and to include a fraction of up to a few percent of the quasar light [see Eqs. (|8|) 
and (^)]. Its contrast relative to the stellar emission by the host galaxy should increase towards 
the rest-frame UV. Current photometric observations of quasar hosts are sensitive to scattered 
light (fuzz) within several arcseconds (^ 10 kpc) from the quasar, at a minimum level of ~ 1% 
of the quasar emission (e.g., McLeod & Rieke 1995; Bahcall et al. 1997, and references therein). 



^While it is possible also to search for broadening of stellar emission lines in starburst or normal galaxies, the 
scattered flux in these cases would typically be much lower than for bright quasars. 
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The method described here requires low-resolution spectroscopy at a comparable contrast level. 
The signal-to-noise ratio of the spectroscopy can be improved considerably by filtering out only 
polarized light from a resolved halo image. Under more favorable contrast ratios, the feasibility of 
detecting extended Lyct emission around high redshift quasars has already been demonstrated in 
many cases (e.g., Heckman et al. 1991a,b; Bremer et al. 1992; Hu, McMahon 8z Egami 1996). 

The scattering in elliptical galaxies might benefit also from the presence of a surrounding 
intra-group or intra-cluster gas. In spirals, the scattering by hot supernova-heated gas would add 
to the more diffuse scattering in the surrounding halo. The effect of hot interstellar gas in galactic 
disks would be heavily contaminated by dust scattering and absorption and by stellar emission. 
Although the scattering by dust should also be polarized, it would not increase significantly the 
width of the quasar emission lines. For example, the broad Mg II emission observed from the 
extended image of the radio galaxy 3C 265 (Dey &: Spinrad 1996) is most likely associated with 
dust scattering of nuclear emission, based on its flux, polarization, and the constancy of the 
measured line width as a function of radius. 

We note that Compton heating or cooling by the quasar radiation itself has only a weak effect 
on the halo gas over a typical quasar lifetime. The characteristic timescale for Compton heating or 
cooling equals ~ 3 x lO^LJ^rfll — 2T'j / U2qq\~^ years, where L46 is the quasar luminosity in units 
of 10^^ erg s~^, and T^ is the Compton temperature of the quasar radiation in units of 10^ K (its 
typical value; see Mathews k. Ferland 1987). This timescale exceeds the Hubble time at z 3 for 
f ^ 2L4^^ kpc. 

The surface brightness distribution of scattered light can be used to infer the density 
distribution of the halo gas. When combined with a measurement of the gas temperature, such 
data could constrain the mass profile of the dark matter as well. This method, however, might be 
compromised by quasar variability (Wise Sz Sarazin 1992). The scattered light from a radius r is 
delayed relative to the central quasar emission by ~ 3 x 10^(r/10 kpc) years, and so variability of 
the quasar on this timescale would affect the brightness distribution of the halo. 

Detection of scattered light can also probe beaming in the quasar emission (Fabian 1989; 
Sarazin & Wise 1993). In the case of beamed broad line emission, the inferred electron temperature 
would be uncertain by a factor of (1 — cos^)^^ due to the unknown three-dimensional orientation 
of the beam relative to the line of sight. Also, the net polarization of the scattered light from the 
entire halo would not vanish in this case as it does in a spherical geometry. Non-sphericity or 
dumpiness in the distribution of infalling gas around the quasar host (Rees 1988) would result in 
similar effects. 

In cases of beamed quasar emission, the contrast between the nuclear and host luminosities 
might be substantially reduced when the quasar beam is not directed at the observer. Indeed, 
the possibility that dust or electron scattering of quasar light accounts for the alignment between 
the radio and optical axes of radio galaxies is supported by observations of some objects (e.g., di 
Serego Alighieri et al. 1989; Dunlop & Peacock 1993; Cimatti et al. 1993, 1994, 1997; Jannuzi et 
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al. 1995; but see also Dey et al. 1997; Pentericci et al. 1998). However, only a small minority of 
all quasars are radio loud and so these objects are by no means representative of the entire quasar 
population. 

The existence of hot gas in halos of high-redshift galaxies without quasars can be probed by 
complementary measurements in the X-ray and radio bands. The forthcoming AXAF^ satellite 
will achieve arcsecond resolution and might detect free-free emission from the brightest halos. 
High-resolution radio telescopes, such as the VLA or MERLIN, could be used to search for the 
associated Sunyaev-Zel'dovich effect of these halos and could provide an independent estimate of 
their temperature when combined with X-ray imaging data. 

Finally, we note that the scattering efficiency of electrons or dust yields a maximum ratio 
between the luminosity of quasars and their hosts. Although the colors and morphology of the 
fuzz around low-redshift quasars indicate association with starlight (e.g., McLeod & Rieke 1995; 
Bahcall et al. 1997), scattering may dominate the fuzz around high redshift quasars at a time when 
galaxies have not converted yet most of their gas into stars. It would therefore be interesting to 
search for an extended scattering component with quasar colors in the photometry of high-redshift 
hosts (e.g., Kukula et al. 1998). The relative significance of starlight versus scattered quasar light 
would indicate whether stars preceded quasars throughout the entire history of galaxies. 

This work was supported in part by the NASA grant NAG5-7039. I thank Daniel Eisenstein, 
Ari Laor, and George Rybicki for useful discussions. 



REFERENCES 

Adelberger, K. L., Steidel, C. C, Giavalisco, M., Disckinson, M., Pettini, M., & Kellogg, M. 1998, 
ApJ, 505, 18 

Bahcah, J. N., Kirhakos, S., Saxe, D. H., k Schneider, D. P. 1997, ApJ, 479, 642 
Bahcah, N. A. & Chokshi, A. 1991, ApJ, 380, L9 

Baugh, C. M., Cole, S., Frenk, C. S., & Lacey, C. G. 1998, ApJ, 498, 504 

Bremer, M. N., Fabian, A. C., Sargent, W. L. W., Steidel, C. C., Boksenberg, A., k Johnstone, R. 
M. 1992, MNRAS, 258, 23p 

Bregman, J. N., & Glassgold, A. E. 1982, ApJ, 263, 564 

CassinelH, J. P., & Hummer, D. G. 1971, MNRAS, 153, 9 

Cimatti, A., di Serego Alighieri, S., Fosbury, R. A. E., Salvati, M., & Taylor, D. 1993, MNRAS, 
264, 421 



http : / /asc . harvard . edu 



-8- 



Cimatti, A., di Serego Alighieri, S., Field, G. B., & Fosbury, R. A. E. 1994, ApJ, 422, 562 
Cimatti, A., Dey, A., van Breugel, W., Hurt, T., & Antonucci, R. 1997, 476, 677 
Dey, A., & Spinrad, H. 1996, ApJ, 459, 133 

Dey, A., van Breugel, W. Vacca, W. D., & Antonucci, R. 1997, ApJ, 490, 698 
Dirac, R A. M. 1925, MNRAS, 86, 825 

di Serego Alighieri, S. Fosbury, R. A. E., Quinn, P. J., & Tadhunter, C. N. 1989, Nature, 341, 307 
Dunlop, J. S., k Peacock, J. A. 1993, MNRAS, 263, 936 
Ellingson, E., Green, R. F., & Ycc, H. K. C. 1991, ApJ, 378, 476 
Fabian, A. 1989, MNRAS, 238, 41p 
Fabbiano, G. 1989, ARA& A, 27, 87 

Fineschi, S., Kohl, J. L., Strachan, L., Gardner, L. D., Romoli, M., Hubcr, M. C. E., & Noci, 
G. 1997, presented at The Corona and Solar Wind Near Minimum Activity, Fifth SOHO 
Workshop, cd. A. Wilson (Noordwijk, The Netherlands: ESA), ESA SP-404 

Fisher, K. B., Bahcah, J. N., Kirhakos, S., & Schneider, D. P. 1996, ApJ, 468, 469 

Haehnelt, M. G., & Rees, M. J. 1993, MNRAS, 263, 168 

Haiman, Z., & Loeb, A. 1998, ApJ, 503, 505 

Hartwick, F. D. A., & Schade, D. 1990, ARA& A, 28, 437 

Heckman, T. M., Lehnert, M. D., van Breugel, Wil, & Miley, G. 1991a, ApJ, 370, 78 
Heckman, T. M., Lehnert, M. D., Miley, G., & van Breugel, Wil 1991b, ApJ, 381, 373 
Hu, E. M., McMahon, R. G., & Egami, E. 1996, ApJ, 459, L53 

Jannuzi, B. T., Elston, R., Schmidt, G. D., Smith, P. S., & Atockman, H. S. 1995, ApJL, 454, 
Llll 

KaufFmann, G., White, S. D. M., & Guiderdoni, B. 1993, MNRAS, 264, 201 

Kukula, M. J., Dunlop, J. S., Hughes, D. H., Miller, L., O'dea, C. P., & Baum, S. A. 1998, 
NICMOS and the VLT: A New Era of High Resolution Near Infrared Imaging and 
Spectroscopy, eds. W. Freudling k R. Hook, Pula, Sardinia, Italy, June 26-27th, 1998, ESO 
Conference and Workshop Proceedings 55, p. 198 

Magorrian, J., Trcmainc, S., Richstonc, D. Bender, R., Bower, G., Dressier, A., Fabcr, S. M., 
Gcbhardt, K., Green, R., Grillmair, C, Kormedy, J., k Lauer, T. 1998, AJ, 115, 2285 

Mathews, W. G., k Ferland, G. J. 1987, ApJ, 323, 456 

McLeod, K. K., k Ricke, G. H. 1995, ApJL, 454, L77 

Mihalas, D. 1978, Stellar Atmospheres, 2nd ed. (San Francisco: W. H. Freeman), p. 420 
Peterson, B. 1997, Active Galactic Nuclei, (Cambridge: Cambridge Univ. Press), §5 



- 9 - 



Pentericci, L., Rottgering, H. J. A., Miley, G. K., McCarthy, P., Spinrad, H., van Breugel, W. J. 
M., & Macchetto, F. 1998, A & A, in press, |astro-ph/9809056 

Rees, M. J. 1988, MNRAS, 231, 91p 

Rybicki, G. B., k Hummer, D. G. 1994, A & A, 290, 553 

Sarazin, C. L., & Wise, M. W. 1993, ApJ, 411, 55 



Somerville, R. S., Primack, J. R., & Faber, S. M., 1998, MNRAS, submitted, |astro-ph/9806228 
Wise, M. W., & Sarazin, C. L. 1992, ApJ, 395, 387 



This preprint was prepared with the AAS lATp^X macros v4.0. 



- 10 - 




Fig. 1. — Profile of a scattered line from a spherical halo [solid line; see Eq. (^)] in comparison with 
the profile for 90° scattering and the same flux normalization (dashed line). Both profiles have a 
unit integral and a unit second-moment. 



